Low-temperature transport of correlated electrons 
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Transport properties of a single-channel Luttinger liquid impinging on a barrier have been studied for g = 1/2 — e, where 
p is the dimensionless interaction constant and |e| -C 1. The relevant diagrams contributing to the conductance are identified 
and evaluated in all orders. Our approach represents a leading-log summation which is valid for sufficiently low temperature 
and small voltage. The asymptotic low-temperature corrections exhibit a turnover from the T 2 ^ 9 ~ 2 behavior to a universal T 2 
law as the voltage is increased. 



PACS numbers: 72.10.-d, 73.40. Gk 



The importance of many-body correlations for trans- 
port phenomena in quasi-one-dimensional (ID) quantum 
wires has attracted a lot of recent interest. Many theo- 
retical studies treat the Coulomb interactions among the 
electrons in the Luttinger liquid framework, and then 
take into account one or a few impurities or barriers 
fi],HJl,§|||. Th ese studies are of direct relevance to 
the tunneling of edge state excitations in the fractional 
quantum Hall (FQH) regime , and to electron trans- 
port through a constriction in narrow high-mobility het- 
erostructure channels Q . Surprisingly few effort has been 
devoted to the nonlinear transport properties of such in- 
teracting ID electron models. In this Letter, we provide 
rigorous results for the nonlinear dc conductance G(V, T) 
in presence of strong interactions, namely near the value 
g = l/2 of the interaction constant. Putting g — 1/2 — e 
with |e| <C 1, we compute the coefficients Aj in the series 
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G(V, T) = > Aj (V, T) e> ln J [max(fc B T, eV)} 



3=0 



(1) 



whereas contributions ~ e 1 ln fc [max(fcBT, eV)] with k < j 
are disregarded. Such a leading-log calculation is valid 
at sufficiently low temperature T and small voltage V, 



e « 1 



eV/k B T K « 1 , T/T K « 1 , (2) 



where Tk is the Kondo temperature. 

If the Coulomb interaction between electrons does not 
contain unscreened long-range tails Q and one is allowed 
to neglect backscattering, the ID electron fluid behaves 
as a Luttinger liquid jl0],[ll| . The interaction strength is 
then fully specified by a parameter g, where g < 1 for 
repulsive interactions and g = 1 for the noninteracting 
case. Introducing a pointlike scatterer to model the bar- 
rier, one finds several remarkable features (l]J^||. Specif- 
ically, for g < 1 the conductance vanishes as T — > 0, im- 
plying that an arbitrarily weak barrier causes complete 
reflection. It is this nonperturbative regime we discuss 
here. Results in other regions of parameter space follow 
from duality relations |12[ or by resorting to perturba- 



tion theory [[[). 

Low-temperature transport is of particular interest 
here. Scaling ideas have been employed by Kane and 
Fisher (KF) to piece together perturbative results in dif- 
ferent parameter regimes. By matching onto the RG 
flows into the stable fixed point, they argued that the 
linear conductance in the nonperturbative regime should 
obey the power law G(0,T) - T" 1 !^ 2 as T —> @. On 
the other hand, there have been given arguments that 
the nonlinear conductance should exhibit a universal (g- 
independent) T 2 enhancement These results being 
correct would imply that the limits V — > and T — > 
do not interchange. The exactly solvable case g = 1/2 
]l|, p^|Jl4| cannot distinguish these two power laws and it 
is essential to study g ^ 1/2. 

Very recently, the KF power law for the linear con- 
ductance appears to have been confirmed experimentally 
0. The tunneling of edge state excitations in the FQH 
regime for filling factor v = 1/3 is well described by the 
transport model considered here with g = v 0, leading 
to a clean experimental setup undisturbed by backscat- 
tering. The findings of Ref. have also received novel 
theoretical support. Integrability arguments can be ex- 
ploited for g = 1/3, and using fermionization || and 
the thermodynamic Bethe ansatz ||, some exact results 
have been obtained. Our study is similar in spirit to 
recent calculations by Matveev, Xue and Glazman ||. 
Their work assumes the presence of weak Coulomb inter- 
actions and therefore corresponds to a leading-log sum- 
mation for g = 1 — e. Unfortunately, it is impossible to 
apply their method around smaller values of g which is 
important because the Kondo scale Tk vanishes with an 
essential singularity as g — > 1~. In contrast, we perform 
a leading- log calculation for g = 1/2 — e. Finally, we 
mention that several quantum Monte Carlo results have 
been given in Refs. [ p^| , ^6[ . 

Within the standard bosonization procedure QJ11]] , the 
fermionic field operator is described in terms of two 
bosonic fields 4>{x) and 9{x) satisfying the equal-time 
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commutator [<f>(x), 6(x')] = (— i/2)sgn(x — x'). Omitting 
irrelevant terms, the Hamiltonian density describing the 
low-energy behavior is ]lj] 



H(x) = |(V0) 2 + l(V0) 2 

+ V 5{x) cos[2^] + eV S(x) . 



(3) 



The (weak) barrier is modelled by a short-range scatter- 
ing potential centered at x = 0, and the coupling con- 
stant Vq is proportional to its Fourier transform at 2fcp. 
The last term represents an external static voltage V. 

There exists a profound connection between this Lut- 
tinger liquid picture and the dissipative quantum diffu- 
sion of a light particle in a periodic potential. We will 
exploit an exact formal mapping between the Luttinger 
model (||) and the dissipative tight-binding (TB) model 
discussed in Refs. jl2];[l3 14 1. The central parameters in 
the dissipative TB model are the hopping matrix element 
A, the dimensionless Ohmic system-bath coupling K and 
the bias a. To avoid spurious divergences, one also needs 
to introduce a cutoff frequency uj c for the bath spectral 
density, which amounts to the bandwidth in the Lut- 
tinger model JhJ. In effect, there is a one-to-one corre- 
spondence between the Ohmic bath modes and the plas- 
mon modes of the Luttinger liquid away from the impu- 
rity. The current-voltage characteristics for transport of 
a Luttinger liquid is thereby expressed in terms of the 
mobility-bias relation for the dissipative TB model. 

This formal equivalence can be shown directly via uni- 
tarily transforming the Hamiltonian @ onto the TB 
Hamiltonian. Employing the standard definition of the 
nonlinear mobility /i(er, T) (iJ|l| and denoting the free 
Brownian mobility by /iOi we find for the current-voltage 
characteristics 

J(V, T)/G V =: G{V, T)/G = 1 - n(Kv, T)/ Mo (4) 

with Go = ge 2 /h. The necessary parameter identifica- 
tions are 
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K 



V = hA 



eV = ha 



(5) 



In the remainder, we discuss the nonlinear mobility for 
K = 1/2 — e (where |e| <§; 1), and then draw conclusions 
concerning the original transport problem from Eq.(Q). 

To compute the mobility, we employ Feynman- Vernon 
theory. One can envision the possible paths of the TB 
particle as moving on a lattice spanned by the forward (q) 
and backward (<?') real-time paths, with the action con- 
taining a nonlocal influence functional due to the elim- 
inated bath modes. Switching to symmetric and anti- 
symmetric combinations of q and q' , the symmetric ones 
can be integrated out, and it suffices to sum over the 
off-diagonal "charges" ^ = ±1 indicating the hopping 
direction of the antisymmetric path. In the end, the mo- 
bility takes the form of a power series in A 2 0,0, 



n(<r t T)/iM) = (2irK/a)ImU(a,T) 
with U describing an interacting kink gas 

°° poo 

U=J2 (-l) m A 2m / dn ■ ■ ■ dr 2m _! ]T 

m=l ^° {£} 

X exp [ - irr^p, „-, + ^ €jS(Tjk)€k\ 



J 



j>k 



2m -1 



x j J sm(Trpj >m K) 

3=1 



(6) 



where the charges £j = ±1 have to obey overall neutral- 
ity, Y) , £j = 0. The interaction potential is 

S(t) = 2Kln[(huj c /Trk B T)smh(7rk B TT/h)} . 

We use the quantities Pj t m = Yji>j & to measure how 
far off-diagonal a given path consisting of 2m charges is. 
With the convention £i = —1, all Pj tTn have to be positive 
integers. Finally, the time interval between the two hops 
£j and is denoted as Tjk ■ 



FIG. 1. Diagram resulting in the contribution to the 
mobility. Time flows along the horizontal line and vertical 
lines symbolize charges £ = ±1. 

For K = 1/2 — e, the phase factors appearing in Eq.(||) 
are 



(_1)(p-1)/2 



p odd, 



7rep 



(_!)(p-2)/2 ; pi 



(7) 



If we consider, e.g., the charge configuration £ = ( 

++), the contribution to U will have a 2ire prefactor. On 
the other hand, the breathing mode integral of the inte- 
rior dipole has a 1/e singularity arising from the l/r 1_2c 
short-time behavior of the intradipole interaction fac- 
tor exp[— S(t)]. Together with the phase prefactor, the 
dipole gives the finite breathing mode contribution 7 with 
7 = 7(w c /7) 2e , where 7 = 7rA 2 /uj c . The frequency scale 
7 defines the Kondo temperature Ik (see below). It is 
then convenient to split up the intradipole time interval 
into a short-time part < 7Y < 1 ("collapsed dipole") 
and a remaining long-time part. 

Next we observe that a collapsed dipole has no interac- 
tions with other charges. Therefore, the grand canonical 
sum over all possible arrangements of collapsed dipoles 
between two confining charges is done easily. Generally, 
there are two types of collapsed dipoles, namely (— +) 
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and (H — ). Because of the phase factors (Q), these can- 
cel each other completely if the exterior charges define 
an even p value (one has to make an excursion to odd 
p which leads to a factor ±1 depending on dipole type). 
However, if one has odd p between the outer charges, the 
respective phase factors for the two types of collapsed 
dipoles are different in magnitude. As a consequence, 
one has to insert a dilute gas of them. With the two ex- 
terior charges spanning a time interval r, this insertion 
simply results in a factor exp(— jt). We shall mark this 
factor by a square in the diagrams below. 



FIG. 2. Diagram giving the contribution U^ 2 \ 



(a) 



grouping the In charges into n pairs where paired charges 
must have opposite sign. By virtue of such a Wick the- 
orem, one obtains n\ graphs per diagram. The sign of 
each graph is determined by the number of crossings if 
one connects paired charges. (4) Group the graphs into 
classes regarding the dependence on even time intervals. 
These are the relevant time integrations which turn out 
to be responsible for leading-log contributions. (5) Eval- 
uate every class in leading-log accuracy [in the sense of 
Eq.(|l])]. The remaining odd time integrations can be re- 
duced to standard integrals. (6) Take into account the 
e-dependence of the interactions again. 



(a) 



r 



(b) 



r 



7=r 



(b) 



FIG. 3. Diagrams resulting in the contribution U^K 

To proceed further, we split up U into terms each 
of which contributes to order e ra_1 (and higher orders). 
The quantity corresponding to the diagram in Fig. |l| 
gives already the exact solution for K = 1/2 Gen- 
erally, the quantity U^ n > is the sum of all diagrams with 
2n charges of total charge zero which cannot be divided 
into neutral clusters. These diagrams are dressed by in- 
serting a gas of collapsed dipoles at every odd time inter- 
val. There are no insertions of collapsed dipoles at even 
time intervals for the reason given above, and these time 
integrations are restricted to 7T > 1 since collapsed con- 
tributions are already included in lower-order diagrams. 
The diagram leading to is shown in Fig. || which 
together with yields all terms in order e. To give 
a final example, one has to evaluate the two diagrams 
shown in Fig. || if one is interested in U^. 

We have extracted the leading logarithms by applying 
the following systematic procedure. (1) Draw all possible 
diagrams contributing to f/( n ). (2) Consider the charge 
interactions S(t) for e = first. (3) Apply a decompo- 
sition theorem to reduce the interactions to pair interac- 
tions. In diagrammatic parlance, this can be achieved by 



FIG. 4. The 2! graphs obtained from the decomposition 
theorem for the contribution shown in Fig. 2. The rela- 
tive signs are + for graph (a) and — for graph (b). The curves 
display the interactions between paired charges. 

To illustrate these rules, the diagram shown in Fig. 2 
can be decomposed into the two graphs depicted in Fig. 4. 
It is then straightforward to find all leading-log contri- 
butions due to these two graphs. Similarly, higher-order 
diagrams can be evaluated in leading-log approximation. 
Computational details will be given elsewhere [p7| . 

It is still a formidable task to evaluate all diagrams 
contributing to the coefficients Aj(a,T) in Eq. ([!]). For- 
tunately, for the asymptotic low-energy properties, it is 
sufficient to calculate them to lowest order in a and T. 
Since Aj(a, T) can be expanded in powers of a 2 and T 2 , 
the leading corrections come from Aj ~ a 2 ,T 2 . The di- 
agrams contributing to these terms can be classified and 
evaluated in all orders in e, yielding the exact asymptotic 
behavior under conditions (0) . The relevant diagrams in 

order n consist of external charge pairs ( ) and (H — h), 

with n — 2 dipoles of finite length inserted in between 
(see Fig. 3 for n = 3). One has again two types of ex- 
tended dipoles and therefore 2 n ~ 2 diagrams in order n. 
We have evaluated the leading logs due to these diagrams 
for any n. All remaining diagrams do not contribute to 
the asymptotic low-energy behavior. 
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For the leading low-energy corrections, we find the ex- 
act result 

H(a, T)/ Mo = 1 - (1/3) [(a/7) 2 + (7Tk B T/hjf] 1+ie ± • • • 

(8) 

If the external bias exceeds the thermal scale, ha ^> k B T, 
Eq.(pJ) gives e-independent low-temperature corrections 
~ T to the nonlinear zero-temperature mobility. On 
the other hand, for ha <C k B T, the asymptotic low- 
temperature corrections are ~ y2+8e_ Thus our dynam- 
ical calculation corroborates the KF scaling law. Eq.(Q) 
with rtq) describes a smooth turnover from the KF law 
G ~ T-/ 9 ~ 2 to the universal T 2 low-temperature behav- 
ior as the voltage is increased. For any nonzero voltage, 
the ultimate low-temperature behavior will be T . 

These findings demonstrate that the limits V — > 
and T — ► do not commute. The anomalous exponent 
1 + 4e = 1/g — I in (||) reflects the critical nature of the 
model at T = V = 0. The turnover from the critical 
7 l 2/g-2 j aw .|- ana iytical T 2 enhancement is due to 
the breakdown of scale invariance when an external volt- 
age is applied fLq] . The physical origin of the T 2 law 
is the low-frequency thermal noise associated with the 
Ohmic nature of the plasmon modes in the leads. Eq.(j|) 
describes the full turnover quantitatively. We wish to 
stress that backscattering (BS) cannot alter the T 2 be- 
havior. In the spinless case, BS is treated as exchange 
event of forward scattering. As this simply leads to a 
redefinition of g, the above results remain valid. In the 
spin-i case, BS gives weak logarithmic temperature cor- 
rections to the T 2 l g ~ 2 law at zero voltage ||. For finite 
voltage, the RG analysis shows that these corrections are 
not present, and one finds again the T 2 law. 

The KF law for the linear conductance holds for tem- 
peratures 

T « T K = (V^e) 9/(1 ~ s) ttFoAb • 

This estimate for the Kondo temperature follows from 
the convergence radius of the series (JsJ) , and for g = 
1/2 — e, the quantity 7 emerges indeed. For T 3> Ik, one 
can safely apply perturbation theory. The T 2 ^ 9 ~ 2 scal- 
ing has been observed in the above-mentioned g = 1/3 
FQH transport experiments . As the source-drain volt- 
age drop was of the order eV jk B T w 0.2, the voltage was 
unable to spoil the KF T 4 behavior. It would be interest- 
ing to perform these experiments at voltage eV 3> k B T, 
where the predicted T 2 enhancement should be observ- 
able. Sample heating is not expected to cause serious 
problems since the conductance remains small. 

Finally, we briefly address the higher powers of a 2 
and T 2 in the coefficients Aj. For the next-leading 
contributions Aj ~ a 4 , T 4 , a 2 T 2 , diagrams of the type 

( X X + +X + +), where X stands for a grand 

canonical gas of finite- length dipoles, have to be consid- 
ered in addition to the diagrams ( X + +) leading 



to Eq.(||). Our results show, e.g., that the next-leading 
corrections in Eq. (|8j) behave as 7 l4+16e (for a = 0). 
In addition, the full asymptotic expansion of the low- 
temperature nonlinear conductance can be deduced in 
closed form |l7| ]. We believe that our leading- log tech- 
nique will also be valuable for exact noise calculations in 
the asymptotic regime. 
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